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The  paper  consists  of  six  sections:  the  first  is  devoted  to  chain 
binomial  methods  and  their  use  in  the  statistical  analysis  of 
measles  and  hepatitis  data.  A  second  considers  time  dependent 
results  for  carrier-borne  epidemics,  and  the  use  of  matrix 
methods  in  computing  probabilities  of  their  final  size.  The 
third  surveys  the  application  of  perturbation  techniques  to 
the  general  stochastic  epidemic,  and  the  estimation  of  infec¬ 
tion  and  removal  parameters  in  this  model  on  the  basis  of  small¬ 
pox  data.  The  fourth  section  sunmarizes  asymptotic  results  for 
the  general  stochastic  epidemic  when  the  initial  populations  of 
susceptibles  and  infoctives  are  both  very  large.  In  the  fifth, 
aorne  recent  results  are  outlined  on  the  costs  of  epidemics;  these 
depend  on  the  stochastic  peth  Integral  under  the  infective  curve. 
Finally,  a  brief  account  is  given  of  the  analysis  of  space-time 
interactions  in  epidemic  processes. 
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1 .  Introduction 

In  hie  basic  paper  on  stochastic  point,  processes  Par'i^l' 
(195*0  first  presented  some  general  methods  for  Mu-kcv  pr  cesses 
( Xt  t ) } j  applicable  to  the  theory  of  epidemics. 

These  methods  consisted  essentially  in  tne  der;  ion  of  a 
symbolic  equation 

M  ’  «I9>  31«>  1  'f  'l-1' 


for  the  characterisoic  function  fid  )  *  £tel9X  •  <  f  process 

X(t),  with  operating  only  on  f.  The  function  ;>  j yf  X)  1.  j 

defined  as 


lira 
5t-T  0 


i0ix(tr&t.  )-X(t  ij 

e _ _ _ 

bt' 


U?i 


gi  ven  that  X(t,>  =  x,  An  analogous  equation  for  the  pro l at  1  iit.y 
generating  function  P(z)  -  ft,-i  in  z J  of  the  pr  >cese  can  also  he 
derived,  which  corresponds  to  the  results  from  the  forward  Kolmogorov 
equation.  The  method  can  be  extended  to  vector  vaiued  Markov  processes 
tX(t)} 

1 
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In  the  context  of  epidemics,  a  point  process  in  time  t  vill 
have  associated  with  it  some  random  vector  X(t)  whose  components 
may  he  one  or  more  of  the  following  random  variables,  susceptibles 
R(t),  infectives  Sit),  carriers  Y(t),  their  locations  (L(t)j  in 
the  plane,  and  other  related  random  variables  such  as  the  cost  Ctt) 
of  the  epidemic.  For  convenience,  epidemic  processes  ore  usually 
aosumed  to  be  Markovian. 

It  is  the  purpose  of  this  review  to  outline  some  of  the  recent 
work  carried  out  in  mathematical  epidemiology.  In  selecting  the 
material  to  be  discussed,  several  valuable  developments  have 
necessarily  had  to  be  omitted.  Readers  interested  in  a  comprehensive 
survey  up  to  1967  should  consult  Dietz  (1967);  a  broad  sketch  of 
current  trends  in  epidemic  theory  is  provided  by  the  17  communications 
in  the  WHO  Symposium  in  Quantitative  Epidemiology  (1971)-  This  is  to 
appear  in  print  very  shortly. 

In  the  present,  paper  I  have,  while  attempting  to  cover  a  wide 
range  of  topics,  been  guided  in  my  choice  largely  by  personal  interests, 
My  hope  1b  that  I  shall  succeed  in  making  those  recent  developments 
which  have  attracted  my  attention  of  aB  much  interest  to  my  readers 
as  myBelf. 

The  paper  consists  of  six  sections:  +he  first  is  devoted  to 
chain  binomial  methods  and  their  use  in  the  statistical  analysis  of 
measles  and  hepatitis  data.  A  second  considers  time  dependent  results 
for  c»rrw-bei-nr  ?p1.dutice.  nnd  the  "se  of  matrix  methods  in  computing 
probabilities  of  their  final  size.  Hie  third  surveys  the  application 


2 


of  perturbation  techniques  to  the  general  stochastic  epidemic,  and 
the  estimation  of  infection  and  removal  parameters  this  model  cm 
the  basis  of  smallpox  data.  The  fourth  section  summarises  asymptotic 

results  for  the  general  stochastic  epidemic  when  the  initial  popula¬ 
tions  of  suaceptibles  and  lnfectives  are  both  very  large  In  the 
fifth,  some  recent  results  are  outlined  on  the  costs  of  epidemics, 
these  depend  on  the  stochastic  path  integral  under  the  infective 
curve.  Finally,  a  brief  account  is  given  of  the  analysis  of  space- 
time  interactions  in  epidemic  processes.  We  now  proceed  to  develop 


these  main  themes. 


2 .  Chain  binomial  methods 


A  simple  advance  in  the  discrete  time  treatment  of  small  scale 
epidemics  has  been  the  reformulation  of  Greenwood  and  Reed-Froat  chain 

binomial  models  as  Markov  chains.  A  brief  outline  of  this  was  first 
presented  by  Gani  (1969)  and  later  developed  in  detail  by  Gani  and 
Jerwood  (1971).  Consider  a  discrete  time  epidemic  process  in  which  the 
latent  period  of  the  infection  13  taken  as  the  unit  of  time*  let  the 
random  variable  denote  the  number  of  infected  individuals  Just 

prior  to  time  t  *  0,  1,  —  ,  who  become  infectious  at  t,  while  Rt 
is  the  remaining  number  of  suaceptibles .  Clearly  R  *  R  +  -3  . 

w  A 

If  0  <  p  =  l-q  <  1  is  the  probability  of  contact  between  any 
2  individuals,  then  in  the  Greenwood  model  where  the  number  of 
infectiveB  at  time  t,  when  these  are  non-zero,  1b  assumed  not  to 
influence  the  probability  of  infection  during  (t,  t+l), 


Pr 


r  !  r  -r  r 

..  In  _  1  t  t  t+1  t  +  1 

[Rt+l  =  t+l'Rt  =  V  =  7~TTr~^''\V.  P  1 

U+-L  t.  U+i. 


(t  B  0,  1,  .  .  .  ) 


(2  1) 


This  clearly  indicates  the  Markovian  nature  of  { R  ] v  For  r  *  k,  say, 

x>  0 

the  tranr^tic, probability  matrix  of  the  process  can  be  written  as 


rt+l 


where  Q  ■  diagll,  q,  ...  ,  q  },  and  P  is  the  remaining  matrix  with 
zeros  in  the  diagonal. 

The  epidemic  process  is  assumed  to  stop  at  time  T  =  t  when 

r  -  r  >  0.  The  probability  of  this  event  Is  given  by 
1  t  "" 


Pr(T  «  t)  =  A^Pt'1QE  , 


(2-5) 


where  A^  -  (0,  0,  ...  ,  0,  1]  is  a  (k+l)  row  vector,  and  E  is  the 
(k+l)  coltrnm  vector  of  unit  elements.  From  (2.3),  the  p.g.f.  of  time  T 
to  termination  of  the  epidemic  is  seen  to  be 


V1  -  ep)'1  0qe 


(0  <  e  <  1)  . 


(2.M 


A  similar  expression  can  readily  be  obtained  for  the  Joint  p.g.f.  of  T 
and  the  final  number  of  infected  cases  at  the  end  of  the  epidemic. 

*n  the  case  of  the  Reed-Frost  model,  the  number  of  infected 
cases  s.  at  time  t  affects  Infection  during  the  interval  (t,t+l),  and 

v 


■  5 


Prlst+i  '  8t+i'  Rt+i  =  rt  ‘  8t.+Jst  ~  V  Rt  = 


rt‘ 

at+i:  (rt  -  W 


n  et  *t+i  9t(v6t+i] 

T-  (1  -  q  )  q 


(2.5) 


From  this,  it  is  clear  that  lSt,Rt)  form  a  bivariate  Markov  chain. 
Techniques  used  in  the  manipulation  of  this  chain  are  similar  to  those 
outlined  earlier,  though  the  matrices  are  now  larger  and  more  compli¬ 
cated. 

The  Greenwood  model  may  be  considered  as  a  Markov  chain  imbedded 
in  a  continuous  time  pure  death  process.  Though  there  are  some  similar¬ 
ities  in  their  characteristics,  the  Reed-Frost  model  does  not,  however, 
correspond  to  the  Markov  chain  imbedded  in  a  simple  stochastic  epidemic. 
But  the  reformulation  of  thece  chain  binomial  models  as  Markov  chains 
enables  us  to  overcome  the  reetrictiveness  of  the  Greenwood  and  Reed- 
Froat  infection  schemes.  For  example,  the  Markov  chair,  imbedded  in  the 
simple  stochastic  epidemic  (for  a  single  latent  period),  or  in  the 
general  stochastic  epidemic  may  serve  equally  well  as  suitable  models. 

We  may  also,  in  the  Greenwood  and  Reed-Frost  type  models,  use  non-homo- 
geneous  Markov  chains  to  simulate  changes  in  the  probability  p  of 
infectious  contact,  possibly  due  to  inoculation.  An  example  of  this  is 
considered  in  detail  in  the  final  section  of  Gani  and  Jerwood’s 
(1971)  paper. 

At  a  more  practical  level,  Bailey  and  Alff-Steinberger  (1970) 
have  used  information  from  basic  chain  binomial  models  to  estimate 
parameters  from  an  associated  continuous  time  Markov  process  with 
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infection  parameter  A.  In  this,  the  instant  of  infection  is  followed 
by  a  latent  period  J,  which  is  normally  distributed  N(p, a  ),  and  is 
succeeded  by  an  infectious  peri  ;d  of  constant  length  a-  On  the  basis 
of  Hope  Simpson's  measles  data,  assuming  either  a  Greenwood  or  Reed- 
Frost  model  in  households  of  two  or  three  people,  Bailey  and  Alff- 
Steinberger  found  that  p,  a  were  of  the  order  of  S  and  7  days 
respectively. 

A  similar  analysis  of  Dr.  K.  Petersonb  data  on  Infectious 
hepatitis  assuming  a  Reed-Frost  model  yielded  values  of  16  and  22  days 
for  p,  a.  It  should  be  pointed  out  that  in  this  caBe,  differences 
in  the  cut-off  point  of  the  data  resulted  in  sizable  differences  in  the 
estimates,  but  a  perfectly  satisfactory  fit  is  obtained  from  the  model 
whether  the  cut-off  point  was  9-10  or  12- 1}  days. 


3 .  Time- dependent  results  for  carrier-borne  epidemics 


An  Interesting  development  In  the  continuous  time  stochastic 
theory  of  epidemics  Is  Gillian  Denton's  (1971)  time-dependent,  solution 
for  the  carrier-borne  Infection  previously  discussed  by  Down ton  (1968). 
Consider  a  closed  population  initially  consisting  of  n  >  l  susceptiblec 
and  a  >  1  carriers  at  time  t  *  0.  If  at  time  t  >  0  there  are 
0  <  r  <  n  susceptibles  and  0<s<n+a-r  carriers,  with  the 
remaining  n  +■  a  -  s  -  r  individuals  removed  from  the  population,  then 
the  transition  probabilities  for  the  process  in  the  interval  (t,  t  +  6t) 
are  given  by 


Pr((r,s)  (r-1,  s+l)) 
k  ttts  6t  +  o(6t) 


when  a  BUBceptible  is  Infected, 
and  becomes  an  unue  tec  ted 
carrier) 


Pr((r, a)  (r-l,e)) 

»  (l-7r)  rs  5t  +  o(6t) 


when  a  susceptible  is  Infected, 
detected  and  removed; 


Pr( (  r, s )  -»  (r,  a- 1) ) 
-  ps  6t  +  o(6t) 


when  a  carrier  is  detected  and 
removed. 


Here  the  infection  rate  is  taken  as  1,  ana  the  relative  carrier  removal 
rate  as  p,  while  the  probability  that  nr.  infective  becomes  an  undetected 
carrier  is  0  <  ir  <  1. 

Hie  probability  generating  function  P(z,w,t)  2  P  (T  )  zr  wB, 

r,  s  rfi 

where  p  (t)  is  the  probability  of  r  susceptibles  and  s  carriere 
rs 


■uaBMkh. 
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at  time  t  >  0,  "iven  that  these  are  respectively  a  and  a  at 
t  =  0,  satisfies  the  second  order  partial  differential  equation 


6P 

at 


=  WlWTT  +  1 


7 .  -  7T  J 


_afp 

alov 


,,  ap 
(1  -  w)p  a^ 


(3-1 


subject  to  the  initial  condition  p(z,w,0)  =  zVwa.  If  -we  write 
n 

P(z,w,t)  =  2  2  **  (''j't)  in  Gani  (1967),  and  take  Laplace  transforms 
r=0 

Fr(w,S)  =  /  e‘  t  fr(w,t)  dt,  Re  9  >  0  , 

with  respect  to  time  t,  (3-1)  can  be  reduced  to  a  set  of  first  order 
partial  differential  equations. 

These  can  be  written  in  matrix  form  as 

A(w)  +  97  =  vaE  ,  (5.2 


where  F1  .  F’(w,e)  =  (Pn(w,fl),  ...  ,  F^v,©)),  E’  =  (1,  0,  ...  ,  0), 
and  A(w)  is  given  by  the  (n+l)  x  (n+l)  matrix 


w(n+p )-p 

- 

-nw(7rw+l-7r) 

w(n-l+p )  -  p 

-(n-1)  w(irw+l-7r) 

w(n-'.+p )  -  0 

, 

I 


„  ifl  n-ra  ->  n+a 
a.0  t  _  „  i-l 


P(n,a)  =  lim -T-  l  n  B  )  (  n  3  }  , 

9-*  0  0  J=0  J  n+1  j»a+l  J  n+1 


which,  with  the  exception  of  the  final  vector  t  is  3.1m  0F(O,d), 

6->0 

then 

a 


P(n,a)  =  {P(n,i)}B  . 


(5.6) 


f 


Further,  the  first,  second,  ...  ,  (n+l)th  columns  of  P(n,a)  correspond 
to  the  vectors  of  probabilities  of  the  total  size  of  the  epidemic  for  a 
initial  carriers  and  n,  n-1,  ...  ,  0  initial  susceptibles  respectively. 

As  is  pointed  out  in  the  paper,  these  results  lend  considerable 
power  to  the  matrix  method  used,  particularly  if  one  ia  concerned  with 
computing  values  of  the  probabilities  of  the  epidemic  size  for  in¬ 
creasing  values  of  the  initial  susceptible  population  up  to  n. 

These  can  be  obtained  straightforwardly  from  a  knowledge  of  the  matrices 
Bj  in  (3.^).  The  result  (3.6)  also  holds,  with  the  obvious  minor 

modifications  in  the  B,,  for  the  case  of  the  general  epidemic 

J 

considered  by  Gani  ( 1967 ) . 
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1+ .  The  general  stochastic  epidemic 


Current  studies  of  the  general  stochastic  epidemic  have  been 
concerned  with  theoretical  developments  involving  the  use  of  perturbation 
techniques,  as  well  as  more  practical  methods  for  the  estimation  of 
parameters  relying  on  electronic  computation. 

An  earlier  application  of  perturbation  techniques  was  made  by 
Bailey  (1968)  to  the  simple  stochastic  epidemic  involving  only 
susceptibles  and  infectives.  His  gave  asymptotically  valid  approxi¬ 
mations  for  a  large  population  of  size  N.  After  a  change  of  time 
variable  T  =  Nt,  the  m.g.f.  M^SjT)  *  of  -the  proportion 

{(T)  of  susceptibles  at  time  T  in  the  deterministic  case  is  known 
to  satisfy  the  partial  differential  equation 


dM, 

a 

5T 


(u.i) 


the  infection  rate  being  taken  as  unity.  Ihe  initial  condition  is 
Md(0,O)  =  where  y  -  (N-n)/n  is  the  initial  ratio  of 

infectives  to  susceptibles. 

In  the  stochastic  model  the  equivalent  equation  for  the  m.g.f . 


is 


Approximating  to  the  first  order  in  N  ,  this  reduces  to 


bH  n/,  ®  s  /  d2M  dM  \ 

Sr- 6(1 -5 e)  > 


(4.2) 
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with  the  same  initial  condition  M(0,O)  =  „  en^/N  aB  before. 

Using  only  the  first  order  perturbation  on  the  corresponding  deterministic 
process,  Bailey  was  able  through  an  eigenfunction  approach  to  obtain  use¬ 
ful  approximations  for  the  mean,  variance  and  epidemic  curve  of  the  simple 
stochastic  epidemic . 

Recently  WelS6  (1971)  has  suggested  an  alternative  perturbation 

method  depending  on  the  moments  v  ■  £(fcr)  of  the  proportion  1  of 

r 

susoeptibles.  These  satisfy  the  differential  difference  equations 


3T  *  r(Vl-vr>  +  ”'i  <-1>J  <3«>  »'J' Vj-l-Vj’ 


(fc.3) 


If  v  I®  expanded  in  powers  of  N”X  as 
r 


v  ■  v 

vr  r 


^  +  v^2V2  +  •••  , 


(4.4) 


and  the  first  two  terns  substituted  in  (^O),  a  simple  set  of  differ¬ 
ential  difference  equations  is  derived  which  can  readily  be  solved. 

An  identical  technique  is  applicable  to  the  general  stochastic  epidemic, 
where  the  population  new  consists  of  susceptibles,  infectives  and 
removals,  l.e.,  Individuals  who  become  immune  or  die.  In  Daniels’ 

(1971)  application  of  perturbation  techniques  to  the  general  stochastic 
epidemic,  a  slightly  different  approach  through  the  cumulant  generating 
function  is  used.  It  is  shown  that  the  first  approximations  to  the 
three  second  order  cuaulaats  of  the  process  satisfy  the  set  of 
differential  equations 


1 

r-*  O 

J _ 

xy 

-y 

-X 

0 

1 - 

H  O 

- - 

d_ 

dt 

3 

-2xy 

+  2 

y 

x-y-p  x 

2K(1) 

‘1 

. 

1 

_ l 

.(xtp)y  _ 

0 

y 

X-PJ 

_ 1 

v'fc-5) 


Here,  x  »  K^,  y  *  are  the  first  order  cumulants,  which  are  the 

solutions  of  the  deterministic  equations  for  the  eusceptlbles  and 
infectives  in  the  process,  p  being  the  relative  removal  rate.  While 
(4.5)  cannot  be  solved  explicitly,  numerical  approximations  are 
obtainable . 

The  validity  of  the  general  stochastic  epidemic  model  has  some¬ 
times  been  called  into  question.  Thus,  an  analysis  of  data  and  estima¬ 
tion  of  parameters  for  such  a  model  would  prove  extremely  valuable.  In 
a  recent  paper,  Bailey  and  Thomas  (1971)  have  made  extensive  use  of  on 
IBM  system  3 60  (Model  4o)  computer  to  analyse  data  from  a  smallpox 
epidemic  of  30  cases  in  a  community  of  120  in  SE  Nigeria.  After  a 
detailed  discussion  of  the  livelihood  functions  based  on  periods  between 
successive  removals,  and  the  number  of  such  removals,  the  ML  estimators 
of  the  infection  and  removal  rates  were  calculated.  These  are 
respectively 


$  .  .00168  +  .00047,  r  =  .162  +  .050  (4.6) 


i 


1 

? 


i 


so  that  the  relative  removal  rate  was  found  to  be  $  -  97  +  22. 

To  cut  down  the  rather  lengthy  time  taken  for  an  exact  determina¬ 
tion  of  the  likelihood  function,  a  gamma  type  approximation  vae  used  for 
the  distribution  of  periode  between  the  (u-l)th  and  uth  removals,  when 
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there  were  8  remaining  infectives,  given  the  paRL  hiatory  of  the 
process.  Computation  using  this  approximate  method  was  then  speeded  up 
by  a  factor  of  18,  without  seriouBly  affecting  the  accuracy  of  the 
estimates.  Further  work  is  to  be  done  to  establish  the  value  of  the 
approximation  for  larger  bodies  of  data. 


5.  Asymptotic  results  for  the  general  epidemic 


In  some  recent  work,  Nagaev  and  Startsev  (1970)  have  obtained 
some  interesting  new  asymptotic  results  for  the  general  stochastic 
epidemic.  Let  n  and  m  respectively  be  the  initial  susceptibles  and 
infectlves  in  such  a  model.  Then,  the  transition  probabilities  of  the 
associated  imbedded  Markov  chain  are  known  to  be 


Pr((r,s)  -*(r+l,s)}  - - HlE -  .  p  , 

p  +  n  -  r  r 

Pr((r,s)  »♦  (r,s+l)l  -  '  p  +°n  .  r  -  > 


(5.1) 


where  r,  s  are  infectlves  and  removals  respectively,  p  being  the 

relative  removal  rate.  The  absorbing  boundary  for  this  Markov  chain 

is  the  straight  line  a  =>  m+r. 

If  v  and  T  are  the  final  size  and  duration  of  the 
n  n 

epidemic  starting  with  n  initial  susceptibles,  it  is  difficult  to 

obtain  explicit  results  for  their  probabilities  directly  from  the  Markov 

chain  formulation.  Clearly,  if  we  define  (X.)  ae  independent  random 

J 

variables,  geometrically  distributed  wJth  probabilities  (p.)  of 
k  J 

BuccesB,  and  8.  *  X.,  then 

.  *  j"  J 


Fr'.vQ  >  k)  -  Pt(Sq  <  m,  8X  <  m+1,  ...  ,  <  m+k) .  (5.2) 


But  though  (5.2)  Is  a  well  understood  probability,  its  explicit  expression 
is  not  easily  obtained.  Nagaev  and  Startsev  have,  however,  succeeded 
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in  deriving  several  asymptotic  results  for  it  when  n  and  m  are 
large;  these  results  depend  on  the  relation  of  m  and  n  as  they 
both  tend  to  infinity.  Typical  of  these  are  the  following: 


Theorem  1.  The  probrhv  ty  of  the  final  size  of  a  general  stochastic 
epidemic  is  given  asym*  .tically,  for  large  n  and  m,  such  that 
m  ■  o(n)  where  a  6  exists  for  which  0  <  6  <  ^  <  1-6  or 
1+6  <  ^  «  o(m),  by 


Prt vn  >  an  *  xbJ  -  ♦(*)  U  +  °U))  » 


(50) 


where 


na 


b 


n 


r  2  5  -|i/2 

K  *%>/<*-!>  ] 

[2  nl/2 

Ml  ♦«j)/(i-o) 


S  (1-°>  -  1 


if  6-  >  1 
n 


if  ^  <  1. 
n 


Further,  in  (5*3)  <2  denotes  the  solution  of  a  +  ^  /n(J-a)  *  0,  and 
®(x)  is  the  standard  normal  distribution  function.. 


Theorem  2.  Let  |p|  «  |m(l  -  where  p  »  o(m)  and  n/m^»  o(p  ^), 

then  in  thia  case  the  probability  of  the  final  size  of  thf  epidemic  is 
given  asymptotically  by 

Prl >  aQ  -  xbj  =  #(x)(l  +  o(l)) 
where  an  «  y§nm,  bn 


I 


(5.4) 


Theorem  5.  Let  >  ft  >  0,  then  the  probability  of  the  final  size 
mn  — 

of  the  epidemic  ie  asymptotically  of  the  Poisson  form 


These  results,  quite  apart  from  their  analytic  interest,  are 
clearly  useful  in  characterizing  the  behavior  of  the  general  stochastic 
epidemic  in  large  closed  populations  with  big  initial  numbers  of  both 
susceptlbles  and  lnfectives. 


6.  The  cost  of  epidemics 

A  different  problem,  initially  broached  by  Encker  (1970)  and 
developed  by  Jerwood  (1970)  in  a  recent  note  and  hia  Ph.D.  thesis 
(1971),  is  the  cost  of  epidemics.  Jervood  specifies  the  cost  of  a  simple 
stochastic  epidemic  in  a  total  population  of  N  individuals  with  i 
initial  infectives,  as 

Ci  *  aWi  +  bTi  ’  (6.1) 


where  a,  b  >  0  are  constants.  In  this,  T.  is  the  duration  of  the 

Ti 

epidemic  and  W.  •»  /  8(t)  dt  the  area  under  the  stochaetic  path 

1  0 

between  0  and  traced  by  the  number  of  infectives  S(t)  at 

time  t.  The  definition  (6.1),  a  linear  combination  of  the  duration  of 
the  epidemic  and  the  total  man-hours  lost  by  infectives,  seemB  reali atic 
though  possibly  too  simple.  Using  the  method  due  to  McNeil  (1970), 
the  Laplace-Stieltjes  transform  of  the  cost  is  found  to  be 


(4.„ 


where  the  infection  rate  is  taken  as  1  for  simplicity. 

For  the  general  stochastic  epidemic,  also  starting  from  i 
Initial  infectives  and  b  susceptibles  in  a  total  population  N,  the 


cost  may  be  written  as 


is 

Cis  *  f  (aS(t)  +  b)  dt  -  aW1#  +  bT^ 


(6.3) 


19 


where  a,  b  >  0  are  again  constants,  S(t)  represents  the  number  of 

infectives  at  tine  t,  and  T^b,  are  respectively  the  duration 

T 

time  of  the  epidemic  and  the  area  /  18  S(t)  dt  under  the  stochastic 

0 

Infective  path  between  0  and 

If  wo  define  the  Joint  Laplace-Stieltja6  transform  <p.  (0  ,9  ) 

16  1  c 

of  a. 

mG  M  •&  T 

<pifl(0l,02)  “  £(e  118  2  iB)  »  *e®i/®2>0>  (6.4) 

then  we  can,  again  using  McNeil's  (1970)  technique,  show  that  these 
satisfy  the  bivariate  difference  equation 

°Wi.  '  0  Vi,.  *  *  •  <«•» 

Here  p  is  the  relative  removal  rate,  and  a.  (0, ,6. )  «  (9  f  0  i_1+8+p) . 

IB  i-  c  12 

We  obtain  formally  that 

(6.6) 

...  ,  N-l) . 

Jerwood  has  examined  in  some  detail  the  distributions  of 

CWis,Tl8},  their  expectations  Mlg  »  g Cwla >,  N1#  “  g(Tlg),  and  hence 

the  expectation  £(C.  )  of  the  cost.  He  has  found  chat  (M.  ,N.  ) 

18  is  is 

satisfy  the  respective  difference  equations 


9. 


n 


-£-  +  I 


s  % 


i 
(  JI 


18  A  ^s  '  p  *ti  7*  % 

(l>  *  1#  •  •  •  $  N-bi  b  ■  0>  1/ 
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(8  +  p)  M1q  =  1  +  P^.1>e  *  aMi+l,B.l 


(s  +  0)  n1d  =  r1  +  +  sNi+i,B-i  • 

The  first  equation  can  be  solved  to  give  results  of  the  form 


(6.7) 


(6.8) 


where  the  p  are  themselves  positive  solutions  of  the  equations 
J 

j ^?1(3)%)  bj  -p  <• -i,  Si . 


Ihe  (M,  }  are  shovn  to  have  the  simple  bounds 
is 

~  <  m 

p  -  is  -  P 


(6-9) 


and  graphical  information  computed  for  the  behavior  of  (M  )  when 

X  s 

1  =  1,  4,  7;  8=3,  has  demonstrated  that  these  bounds  are  reasonably 
good. 

The  equation  for  the  {N^}  yields 


(6.io) 


where  N^Q  *  2  k  /p  •  Hence  they  can  be  evaluated  recursively  without 


kd 


difficulty. 

A  large  number  of  interesting  problems  remain  unsolved  in  this 
area;  several  of  these  are  currently  under  investigation. 
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7 .  Space-time  Interactions  In  epidemics 

A  problem  of  great  importance  in  epidemiology  has  been  the 
identification  of  clustering  among  reported  cases  of  a  disease.  Knox 
(1964  a,b)  investigated  methods  for  deciding  whether  leukaemia  cases 
occurring  at  small  distance  from  each  other  also  occurred  close  in 
time;  such  space- time  clustering  would  support  the  theory  that  the 
disease  was  contagious.  His  method  was  later  improved  upon  by  Barton 
and  David  (1966);  a  further  extension  by  Pike  and  Smith  (1968)  allowed 
for  limited  periods  of  susceptibility  and  infectiousness. 

Barton  and  David  specified  Knox's  procedure  in  graph- theoretic 
terms.  Let  X  be  the  number  of  pairs  of  caBes  from  among  the  *, 
reported  which  are  found  to  be 

(i)  living  within  6  fixed  distance  B  from  each  other 
(,ii)  infected  within  a  fixed  time  period  t  of  each  other. 
Adjacency  matrices  S  and  T  of  size  n  x  n  may  be  written  such  that 


S .  »  1  if  case  J  is  at  a  distance  no  greater  than  8  from  case  i 

*  V 

and  is  0  otherwise,  wnile  T^j  =1  if  case  J  occurs  within  a  period  t 


of  .1  and  is  0  otherwise.  By  definition,  8  ■  T  »  0.  The  measure 

J  J  J  J 

X  of  space-time  clustering  is  then  given  by 


(7-1) 


This  integer  is  found  to  be  distributed  approximately  as  a  Poisson 
variable  and  can  consequently  be  used  as  a  test  statistic.  Details  of 
other  gjaph-theoretlc  structures  of  use  in  epidemiology  way  be  found  in 
Tautu  ( 1970’ . 
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In  Knox's  method  the  space  and  time  measures  related  to  cases  i 
and  J  are  either  0  or  1.  But  S. ,,  T .  .  need  not  necessarily  be  re- 

XJ  X  J 

strioted  to  these  values;  any  reasonable  measure  of  contact  in  space 
and  time  may  be  used.  The  variable  X  will  the.i  remain  an  indicator 
of  space-time  clustering;  its  first  two  moments  in  this  general  case 
have  been  obtained  by  Barton  and  David  (1966).  Mantel  ( 1967 )  later 
developed  a  generalized  regression  approach  to  the  clustering  problem 
of  which  the  Knox  and  Barton  and  David  results  are  special  cases. 

Here  a  statistic 

(7.2) 

is  defined  where  S.  R,  are  respectively  functions  of  the  space  and 

X  J  X  J 

time  coordinates  of  cases  1,  J.  Testing  of  the  statistic  is  curried  out 
by  Monte  Carlo  methods,  or  in  some  cases  where  it  appears  Justifiable, 
by  the  use  of  an  approximate  normal  distribution  for  U  *  [Z-  5(Z))/o^. 

In  a  recent  paper  Klauber  (1971)  considers  an  extension  of 
Mantel's  method  to  test  clustering  between  two  sets  of  points.  His 
methods  are  to  be  applied  to  the  diagnosis  of  lymphatic  leukaemia  cases 

in  both  man  and  pet  cats.  Let  (X^,  Y^,  T^)  (i  1,  2,  ...  ,  I)  and 

. 

(X,,  Y  ,  T  J  (J  -  1,  2,  ...  ,  J)  represent  the  space  and  time 
J  J  J 

cooidinates  of  two  sets  of  observations,  say  on  men  and  cats  respectively. 
Ihen  we  may  define  the  space  and  time  distances  between  case  i  among 
men  and  case  J  among  cats  by 
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SiJ  =  S([(xi  “  XJ}  +  (Yi  "  Yj)2]1/2)  »  \2  -  R(|T4  -  T*|)  (7.3) 

where  S,  R  are  given  functions,  for  example  the  reciprocals  of  the 
distances  between  cases,  plus  some  adjustable  constant. 

The  statistic  (7.2)  may  still  be  used*  it  can  be  compared  to 
its  randomization  distribution  by  assuming  one  set  of  points  fixed, 
and  the  coordinates  of  the  other  randomly  permuted,  or  both  sets  may 
be  taken  as  random.  The  mean  and  variance  of  Z  can  be  found  for  both 


randomization  models,  and  in  some  instances  the  approximate  standard 
normal  variable  U  =  (Z  -  §(Z))/cz  may  again  be  used  as  the  test 
statistics.  As  might  be  expected,  this  approximation  improves  for 
larger  samples.  Two  examples  are  given  of  tests  based  on  empirical 
data;  one  of  them  is  concerned  with  lymphoma  cases  for  117  catB  and 
93  dogs  in  California,  where  no  statistical  significance  at  a  reasonable 
level  wbb  found  for  clustering. 

It  will  be  apparent  from  this  brief  sketch  of  current  work  in 
epidemiology  that  the  field  is  in  a  state  of  rapid  development.  The 
intrinsic  interest  of  the  problems  considered,  and  the  liveliness  of 
research  in  this  area  will,  I  hope,  encourage  more  probabilistB  and 
statisticians  to  contribute  to  it  in  future. 
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